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Abstract 

Magnetic fields in the range < eB < 1 GeV 2 effectively probe 
internal quark structure of chiral mesons and test basic parameters 
of chiral theory, such as (qq),fw We extend previously formulated 
chiral theory, including quark degrees of freedom, to the case of strong 
magnetic fields and show that the quark condensate \{qq)\ u ,d grows 
quadratically with eB for eB < 0.2 GeV 2 and linearly for higher field 
values. These results agree quantitatively with recent lattice data and 
differ from \P T predictions. 



1 Introduction 

Strong magnetic fields (MF) are expected to be present in the early universe 
PQH1E], in neutron stars (magnetars) [HE], and also in non-central heavy-ion 
collisions [6]. It is very interesting how MF can modify strong interactions, 
deduced from QCD. In particular, it was recently shown [7J, how MF influ- 
ence asymptotic freedom behavior of the QCD strong coupling constant and 
gluon exchange interaction. The interplay of confinement and MF in neutral 
mesons was studied in [8], showing a dramatic change of masses with growing 
MF. 

Of a particular interest is the influence of MF on the chiral symmetry 
breaking (CSB) dynamics, and in a more general setting, on symmetry break- 
ing pattern in field theory. 

It was emphasized in [9], that MF has the tendency to strengthen the 
chiral condensate, the phenomenon called there the magnetic catalysis, see 
[TO] for a review. 
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The behavior of chiral condensate in MF was studied in the chiral pertur- 
bation theory (xPT) [III [121 EH] , in Polyakov-Nambu-Jona-Lasinio (PNJL) 
[T5] and in other models [16]. On the lattice the corresponding analysis was 
done for quenched QCD in PU Q3], in n f = 2 QCD in [191 [20] and in n/ = 4 
SU(2) theory in [21]. In all cases both in models and on the lattice the chiral 
condensate grows with MF, however in different way. 

Recently a comprehensive analysis of this problem on the lattice with 
physical quark masses was performed in [22] and has shown a quadratic 
growth of condensate up to eB ~ 0.2 GeV 2 and approximately linear behavior 
above this value. 

This linear behavior contradicts the (xPT) results, see e.g. [TTj[T4]. where 
linear growth was found with the slope almost twice as small as in [22], and 
also contradicts PNJL quadratic growth at large eB |15j . 

As it was mentioned in [13], the reasonable region of xPT application is 
for eB < m%, where the quark condensate grows quadratically, in qualitative 
agreement with lattice data however the linear behavior up to 1 GeV 2 is 
outside of xPT reliable region. 

It is the purpose of the present paper to study quark condensate in MF 
starting from the basic QCD equations, derived earlier without MF in [23j 
[2U [221 [2S]. It was shown there, that GMOR relations and expressions for 
(qq) and f w can be derived from the basic QCD quantities: string tension, 
a s and current quark masses in good agreement with experiment and lattice 
data. 

Recently [21] these results were extended to account for growing cur- 
rent quark masses m q and in particular the dependence of (qq) on m q was 
established to be in agreement with lattice data [28] . 

In the present paper we follow the same line of the formalism of [231 [2U 
[251 12H], but now in MF and find the behavior of (qq) for u and d quarks and 
their average with growing MF. 

As a result we observe in the resulting dependence of the average (qq) 
and (uu), (dd), the same features and good quantitative agreement with the 
lattice data obtained in [22]- The physical reason for this dependence of (qq) 
on eB is clarified below in the paper. 

The paper is organized as follows. In the next section a general derivation 
of GMOR relation and expressions for (qq) and /„- are given, in section 3 the 
MF dependence of basic terms is established, in section 4 results are discusse 
and prospectives are given. 
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2 Derivation of the effective chiral lagrangian 
in MF 



We follow here [21] to write first the effective lagrangian of the light quark in 
the confining field of the antiquark in MF, starting with the standard QCD 
partition function in Euclidean space-time 

Z = J DAD^D^ + exp[L + L x + L int ] (1) 

where 

L ° = ~\J **W> ( 2 ) 



Li = -% I f tp + (x)(D + m f ) f ip(x)d 4 x (3) 
L int = I f ^ + (x)gA(x) f tf;(x)d 4 x (4) 



and 

D = ^(d,-ieA^(x)), AW(x) = i[xxB]. (5) 

Averaging Z over vacuum gluonic field and keeping only lowest (bilocal) 
correlators of color fields D^ uM (x,y) = ■^tr(F liU (x)^(x,y)Fx tr (y)^{y,x)), 
one finds as in [241. 



(Z) A = J Lhl>D$ + exp(L 1 + L%), (6) 

where 

Leff = J d A xd A y s ^ a (x) f i) bp (x) 9 i) bl (y) 9 i) aS (y)J aPnS (x,y) (7) 

and 

Jafinsipiy) = {l, J ) a f}{lv)~ 1 5J l ,v{x,y) (8) 



Jy, v {x,y) = 9 J c du a dvpD a ^ v (u,v) (9) 

Here indices f,g refer to flavor, a, b to color and a, (3,[i, v to Lorentz 
indices. Eq. (J5J) implies that some contour gauge is used for simplicity, 
but the final result is gauge invariant and the most important property of 
Jfj, v {x, y) is that it is proportional to the distance of the average point f^rrH 



to the contour C (linear confinement), and the effective distance between x 
and y (nonlocality) is of the order of the vacuum correlation length A ~ 0.1 
fm. 

In the large N c limit the four-quark expression in i^// can be replaced 
by the quadratic one, using the limit 

f ^(x) 9 ^(y) ->■ 5 fg N c f S Pj (x, y), (10) 

where ^S^(x,y) is the quark propagator. 
As a result one obtains the form 

Lfff -i J d A xd A y f^ a (x)^M aP (x, y) ^(y), (11) 

f9 M aS (x,y) = -U^(x,y)(^ f9 S{x,y) lu ) aS , (12) 
and the quark propagator satisfies the equation 

{-ib - im f ) f S(x, y)-if if9) M{x, z) 9 S{z, y) = 5 {4) {x, y). (13) 

It is convenient to use the following parametrization of *M(x,y) in terms of 
scalar functions, flavor singlet M s (x, y) and flavor triplet 4> a (x, y),a = 1, 2, 3, 

M { J 9) x, y) = M s (x, y) exp(i 75 ^ a (x, v)tf = M s(*> v)U$\x, y) (14) 
As a result the effective Lagrangian assumes the form 

L+ = J d^d*y{ty+(x)[(i^ 9 ^{y)} , 

(15) 

and the partition function can be written as 

{Z) A = J DipDiP+DMsDfaexpLf. (16) 

Integrating over DtpD^ + one obtains the effective chiral Lagrangian (ECL) 
Lecl, 

(Z) A = J DM s D(f) a exp L EC l, (17) 

where 

L ECL = N c triog[(iD + im f )l + iM s U\. (18) 
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Finally,the ECL at the stationary point in the integral (I17p is defined by 
conditions = = 0, which yields 

zMf {x,y) = (^S-yJJ^y), # = 0, (19) 

and 

S = S+tfa = 0), 5 = -[(if) + im f )l + iM s U}. 
Insertion of f lT9|) in f fl8|) yields the effective action for pseudoscalar fields 



Lecl -»■ -W{<j>) = N c tr\og[(iD + im f )l + iM s U}. (20) 

Our final step here is the local limit of J fMU (x,y) and Mj®(x,y) proved in 
[23] . which yields 

<f> a {x,y)^<f> a {x), Mf\x,y)^M^){x)5^{x-y) (21) 

Expanding W(4>) in powers of 4> a and keeping quadratic terms, one has 

1 r r1 4 kr1 4 k' 

^ (2) (0) = ~ / ^Ly Uk)N{K k ' )Mn (22) 

where 

AT r 

N(k,k') = dxe^ k+k >tr(AM s ) xx + 



i ^ k+ 2 k \ x + y) + ^{k- k')(y - x) 



J d 4 (x - y)d 4 f— y^) ex P 

x tr[A(x,y)M s (y)A(y,x)M s (x)}, (23) 



with the definitions 

A= (D + m + M s )~\ A=(D-m-M s )~ l . (24) 
Using translation invariance of traces in ( l23l) one can rewrite it as 



N r - r] 4 k 

W^(<P) = fj Uk)U-k)N{k)^ (25) 
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M(0) 

Figure 1: Quadratic terms in the pionic action (18) 



where 

N(k) = hr{{AM s ) + f d 4 ze ikz A{0,z)M s {0)A{z,0)M s {0)} = hr{AM). 

(26) 

The graphical representation of W^ 2 \(f>) is given in Fig.l and the picture, 
which helps to understand the physical meaning of M(x, y) is given in Fig. 2. 
As a result the standard chiral condensate in the Minkowskian space-time 
can be written as 

(#>m = i(ip^ + ) E = -N c trA = (uu + dd) M . (27) 
This can be rewritten as 

trA = -(M(0) + m)G {0 \k = 0), (28) 
where m = m »+ Wd ; and M(0) is expressed via the vacuum correlation length 

A ED 

M(0) = -^=<t\{1 + 0(a\ 2 ) (29) 
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Here A is the vacuum correlator length, defined from the mass in the field 
correlator [29], 0.1 f m < A < 0.2 fm. Therefore 100 MeV < M(0) < 200 
MeV. 

Writing the Green's function G^(k) as a spectral decomposition, one 
obtains 

G(°)(A;) = -£-|U, (30) 

In a similar way, as shown in [23] and [30], the qq Green's function 
G^ MM \k) with M(0) in both vertices is connected to decay constant f n , 
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q y M^x , y ) 




q 

Figure 2: The flux-tube operator M(x, y) in the Wilson loop 



namely 

q(mm)ma = _ M 2 (0) y c l G {MM) {k) - G (MM) (0) = — + 0(k 4 ). 

„ k 2 + 2n c 

(31) 

One can express c n as c n = ^/^■V'n(O) an d finally one obtains chiral relations 
for 5 = 0, see [21] for details of derivation 

N ib 2 (0) 

-(gg>, = AT c (M(0)+mO£^^e- m " A , % = u, d, s (32) 
/* = iV c M 2 (0) £ ^ie" m " A (l + m n A). (33) 

n=0 m ™ 

Note, that the GMOR relation [31J is modified, 

22 mM(0) - _ m„ + m d 

m ' / ' = M(0) + m l(MM> + WI> m = ^~ ' (34) 

however for m — > 0, (jM]) coincides with the original form. 

Eqs. (f3"2"j) . (f3"3"|) express quark condensate and /„. in terms of M(0) and 
the reference spectrum in the PS channel, i.e. the spectrum, where chiral 
symmetry is not taken into account, but only confinement term M s (x) is 
present, the corresponding masses m n , calculated in [24J are m = 0.4 GeV, 
mi = 1.35 GeV, m 2 = 1.85 GeV. Taking A = 0.1 fm and M(0) = 0.15 GeV, 
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one obtains (cf appendix of the second paper in [21]) 

-M = (217MeV) 3M(0>+ ™, /, =9 6GeV M( ° )+r " (35) 
n f K ' (150MeV)' J (150MeV) V ' 

One can see, that our values are in the correct ballpark, close to average 
lattice and phenomenological values. 

Now we turn to the case of nonzero MF. Eq. (1321) in case of MF can be 
rewritten for each flavor separately, 

\w)i\ = N c (M{Q)+rrH)2^\ A 7^5 — + - ttT) — (36) 

n =o \ m n>i m nji J 

where i — u,d,s and the subscripts (H — ) and ( — h) refer to the quark and 
antiquark spin projections on the MF B, and the coefficients | are due to 75 
in the vertices of the Green's function G^°\k), (75 1 — > ^ + ~^ h ^ . 

The main problem reduces to the calculation of the spectrum of eigenval- 
ues and eigenf unctions ip n ,m n , which are to be found from the Hamiltonian 
H containing MF B, and derived from the path integral Hamiltonian [S] 

Ho = U\ — m — + ^ — ^rJ (37) 

which can be rewritten in the form (i = u, d, s) 

^h{-i-i^))^-^ (38) 

where u = JaupL.. 

The eigenvalues m n ,i of Hi depend on tJi,U2, and the final eigenvalues 
m n i entering in (J3J) are obtained as stationary points in variation over ui,U2 



dm n4 



= 0, e.g.mi7 ) (4 0) ,4° ) )^mg-). (39) 

,(0) 
-1.2 



To avoid purely numerical calculations, one can simplify the hamiltonian 
fl38|) . replacing linear confinement V COIt f(r)) in fl38l) by a suitable quadratic 
form, with coefficient found from the stationary point condition 



rrUin) x T r(quadr) & ( 2 M ^ m n,i 



= (40) 

7=70 
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To check the accuracy of this replacement, one can compare the lowest 
mass eigenvalues and i/j n (0) first without MF. In particular, 

while mjj n = f^fp") ^ , and these two figures differ by 4.5%. 
For ipn{0) the corresponding results are 

I^(0)IL=^^ = 0.065a 3 / 2 ; 



l^o(0)| 2 uadr = = 0.065a 3 / 2 (c = 1) (41) 



0-3/2 ("a) 3/4 

This coincidence of |"0o(O)| 2 will be of special importance in what follows, 
since the main effect of MF, as will be seen, is the increase of |^ n (0)| 2 in 
dSH]) due to MF. It is clear, that with growing B the size of meson is de- 
creasing, and the difference between and V^jr will be even more 
suppressed, since these both interactions vanish simultaneously, v}^(r — > 
0) = K ( o q nf adr) (r^0)=0. 

For us it will be most important how |^„(0)| 2 depends on MF, and es- 
pecially, how MF enters in the expansion of |^ n (0) | 2 in powers of B. To 
this end one can see in the Hamiltonian fl38|) . that MF enters via the term 

Vb = (^p-^j ■ It is clear, that in the perturbative series expansion the 
MF enters as (eiB) 2 

M0)b = Vv(0) + 0(( ei B) 2 ) + ... (42) 

As a result of (14*0]) one can immediately write the analytic expressions for 
ip n (0) and m n of the following form, e.g. 



/ 2 \ - 1 / 4 

-i/2 ( a z c^ 



W")r = ^. ri = 2[(e q BY + a 2 c) , r 8 = ^— j (43) 

where c = — , and 

_ m 2 + a; 2 - e.Bo-! m 2 + u; 2 + e g B(T 2 



' 1 ^-B 2 + a 2 c(2n ± + 1) + Va 2 c(n 3 + J)] + (45) 



where 



Expressions ( )36|) . ( )43|) . ( 1441) . ( |45|) contain all information necessary to com 
puter quark condensate for varying MF. 
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3 The MF dependence of quark condensates 

We are interested in eigenfunctions and eigenvalues of (+— ) and ( — h) states 
both for u and d quarks. For the (H — ) state, one can see in (|45|) . (jHJ) that 
at large B, m|,^ 3 ' tends to a constant limit, together with u\ = 002- In this 



case the parameter c in (T43]) and ( j45J) which can be expressed in general as 

c + - = ^= Q) 4/3 4//3; >y{B) = P{B){au,)-W uj{B) = a{B)^ (46) 

and P{B), a(B) are changing in finite limits and one obtains that c H (B) « 1, 

for all B, and m H tends to a constant limit for B — > 00. Hence one can 

write for lowest levels n 3 = 0, 1, 2, ... 



vWe 2 ,^ 2 + a 2 



ivfcUo)i 2 - • v (2 ; )3/2 • ( 47 ) 

Note, that in the limit B — > this expression for |?/>(0)| 2 yields equal values 
for 713 = 0, 1, 2 as it should be for pure linear confining interaction. 

For the ( — h) case the situation is different, and at large B the stationary 
point value Uq + ^ ~ J2\e q \B+ + ^, and the parameter c_ + is increasing with 

c_ +(B )^(l + ^) 2/3 . (48) 



As a result, the |^ (0) | 2 can be written as 



l^; +) (o)r = (- 2 c- + ) 3 / 4 



\ 



e q B\ 2 1 



<j J c__ 



(49) 



Moreover, mt^ « 2,/2|e aJ B| + f, and K c Tt )3/4 = at large 5. 
Combining fj45|) and fj4*9|) we therefore can write 



where c_ + is given in ( )48|) . 
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Table 1: Values of A — u,d given by (15171) in comparison with lattice 

data from [22] 



eB 





0.2 


0.4 


0.6 


0.8 


1 


(GeV 2 ) 



















0.156 


0.48 


0.865 


1.273 


1.65 


this paper 














AE„ [21] 





0.185 


0.51 


0.86 


1.235 


1.60 


lattice 



















0.048 


0.158 


0.308 


0.48 


0.67 


this paper 














A E d [21] 





0.095 


0.23 


0.40 


0.57 


0.73 


lattice 















We shall be using notations of [22J for the increment of quark condensate 
function of MF (denoted n(B) in [9]) 

The resulting values of A^(B) from (jSTIj) . (15T1) are compared in Table 1 
with the corresponding lattice calculations in [22]. 

One can see a reasonable agreement between our theory and lattice data. 
One can simplify the B dependence of Eq. (150]) . writing 



(52) 



Kw)i(B)l = Kw>i(o)|. 



The form f[52|) also satisfactorily describes data of [22] and [20], taking 
Mj as a fitting parameter. In the case of [20]the fitted values of Mj are 
approximately 

Ml « 0.35GeV 2 , M d 2 w 0.27GeV 2 , (53) 

i.e. are 2 and 1.5 times larger, than a (the corresponding c = c + _ = 
c_ + = 4 and 4.25 for u and d cases), however the agreement with lattice data 
is worse. Qualitatively the same situation (with even larger M 2 )takes place 
in comparison with [T7] [T8] . 
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As a whole the behavior (I50j) . (l52p correctly reproduces the main quali- 
tative features of the quark condensate as a function of eB: the quadratic 
behavior proportional to (eiB) 2 at small B and linear behavior ~ \eiB\ at 
large B. It differs from the results of others approaches. In particular, xPT 
[TT] predicts linear behavior in eB with the slope much smaller, than in 
lattice data [22J. 

As it was mentioned in [20J , for A J2 U an d A J2d there is a simple relation 

A E(!)=AE(£), (54) 

u V Z J d 

which is satisfied in lattice data, and, of course, is trivially satisfied in our 
definitions. 

Writing for large B ^> a, that A = aiB, one immediately obtains 

from fl54|) . that a u = lad- This relation is approximately satisfied in lattice 
data [201 [22], an d in our expression fl50|) . fl52l . 



4 Discussion and conclusions 



We have used our formalism for chiral dynamics, presented in [231 [2 
26[[27], which is derived not from purely symmetry considerations, but from 
the QCD quark dynamics, where chiral symmetry appears approximately in 
the small m q limit of the effective QCD Lagrangian. In this way all basic 
degrees of freedom are connected to the confined quarks, and it is finally the 
confinement, which dictates properties of fundamental chiral quantities (qq), 
f n etc., and gives them numerical values, expressed via a. This is in contrast 
with standard chiral ideas, where chiral constants are postulated from the 
beginning and chiral symmetry breaking is not connected to the basic QCD 
quark dynamics. 

Actually our reasoning for the calculation of quark condensate is very 
simple. After one derives scalar confining interaction M s (x), which acts on 
each quark or antiquark, one can write quark condensate (qq) = —N c trA as 



trA — ( tr- 



1 



M s + m q + dl 
(M s (0) + m q ) Ur 



tr 



t:(M s + m„- d) — 

M s + m a + d (M s + m q -d)_ 



1 



75 



1 



^75 — 

M s + m a + d M s + m + d 



(55) 
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since the term with d is odd and vanishes. 

The last quantity on the r.h.s. of (155]) is the qq Green's function, propor- 
tional to ^(0) for each n state. The only difference in presence of MF is the 
replacement d — > (d — ieA^), which immediately gives the proportionality 
(qq) ~ tp 2 (r = 0; eB). This latter quantity is linearly rising with \eB\, since 
MF is "focussing" quark- ant iquark system at small distances. 

External magnetic field is here crucial for our understanding of chiral 
dynamics, and using MF one may decide, what is the role of quark dynamics 
in the chiral phenomena. In this respect the comparison of the behavior of 
quark condensate (or AJ2{B)) as a function of B in different models and 
lattice data is showing the following: 

1. Standard chiral theory at large B ~ o gives linear behavior (qualita- 
tively correct) but with the wrong slope. 

2. Our approach, Eq. ( 13 6 p predicts quadratic behavior at small B B < a 
and linear at larger B with slopes different for u and d quarks, both in 
agreement with existing lattice data. 

3. The PNJL model as shown in [22], is also in disagreement with lattice 
data at larger eB. 

Physically, it is clear, that MF, acting on quark charges, immediately 
discloses the internal quark structure of PS mesons, while the effective chi- 
ral Lagrangian, as in [111 112] describes only internal multipionic degrees of 
freedom. The latter can be important only for B <C l/ r o> where r ~ 0.6 fm 
is the pionic radius, i.e. for B <C 0.1 GeV 2 , while for large B the standard 
chiral picture is irrelevant, as it is confirmed by lattice data. 

These considerations suggest the idea, that the true chiral dynamics can 
be derived e.g. from the effective Lagrangian (JTHjl . and should finally display 
coexisting quark and chiral symmetric degrees of freedom, demonstrating 
how the latter disappear (suppressed) for growing quark masses. How it 
happens with the pionic spectra, was explained in [27], demonstrating the 
unifying spectrum based on chiral and quark degrees of freedom at the same 
time. 

The author is grateful for discussions to N.O.Agasian, M.A.Andreichikov 
and B.O.Kerbikov. 
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